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Abstract — This  paper  has  derived  a  new  method  for  estimat¬ 
ing  the  maximum  ventricular  elastance  Emax  by  assuming 
that  the  following  two  assumptions  are  simultaneously  es¬ 
tablished  in  a  certain  part  of  the  ejection  period:  the  ven¬ 
tricular  elastance  E(t)  changes  linearly  with  time  and  the 
dead  volume  Vq  does  not  change.  The  method  can  estimate 
Emax  in  a  single  beat  without  any  change  in  the  cardiac  load 
by  measuring  the  ventricular  pressure  and  either  the  ven¬ 
tricular  volume  or  outflow.  The  method  does  not  need  to 
previously  assign  the  approximation  range  of  E(t)  and  does 
not  include  the  recursive  procedure  to  search  for  the  peak 
time  tmax  that  gives  Emax •  In  an  experiment  in  vivo ,  it  could 
be  ascertained  that  the  method  can  derive  the  estimate  of 
Emax  roughly  close  to  the  true  value. 

Keywords  —  Emax>  estimation,  end-systolic  ventricular  elas¬ 
tance,  ESPVR,  cardiac  function 

I.  Introduction 

The  maximum  value  (Emax)  of  the  ventricular  elastance 
( E(t ))  has  been  thought  to  be  a  good  index  for  quantify¬ 
ing  the  ventricular  contractile  performance  independently 
of  the  cardiac  load.  However,  the  conventional  method  for 
estimating  Emax  requires  a  change  in  the  cardiac  load  dur¬ 
ing  several  beats.  To  promote  clinical  application  of  Emax, 
it  is  desirable  to  estimate  Emax  low-invasively  on  the  basis 
of  information  acquired  within  a  single  heartbeat  without 
any  change  in  the  cardiac  load. 

As  one  of  such  methods,  the  authors  [1],  [2]  have  pro¬ 
posed  the  parameter  optimization  method  (POM).  The 
POM  can  estimate  Emax  on  the  basis  of  ventricular  pres¬ 
sure  and  either  ventricular  volume  or  outflow  acquired  only 
in  a  single  beat.  The  POM  does  not  need  the  normalized 
time-varying  ventricular  elastance  curve  (-E)v(Uv))  used  in 
Senzaki’s  method  [3].  In  the  POM,  however,  it  takes  much 
computational  time  to  search  for  the  optimal  parameters 
determining  the  linear  function  to  approximate  the  elas¬ 
tance  E(t),  and  the  approximation  range  must  be  given 
previously. 

To  overcome  these  defects  of  the  POM,  in  this  study,  a 
new  method  for  single-beat  estimation  of  Emax  has  been 
developed  by  directly  using  linearity  of  ventricular  elas¬ 
tance  in  the  ejection  period.  In  the  proposed  method,  it  is 
expected  that  the  computational  time  can  be  reduced  and 
that  the  problem  of  how  to  determine  the  approximation 
range  of  the  elastance  curve  can  be  avoided. 


II.  Methods 

A.  Conventional  Method 

Let  P(t)  and  V(t)  denote  ventricular  pressure  and  vol¬ 
ume,  respectively.  Ventricular  contractility  can  be  ex¬ 
pressed  by  the  ventricular  elastance  E(t)  calculated  from 


E(t) 


Pit ) 

V(t)  -  Vo 


(1) 


where  Vo  is  the  so-called  dead  volume  defined  as  the 
volume  axis  intercept  of  the  regression  line  of  the  end- 
systolic  pressure- volume  relation  (ESPVR).  Drastic  change 
in  preload  or  afterload  is  necessary  to  obtain  sufficiently 
distinct  end-systolic  points  to  determine  the  ESPVR.  Emax 
is  equal  to  the  maximum  value  of  E{t). 


B.  Assumptions 

In  this  study,  the  following  two  assumptions  are  pro¬ 
posed. 

Assumption  I:  E(t)  changes  linearly  with  time  in  a  cer¬ 
tain  part  of  the  ejection  period. 

Assumption  II:  Vq  in  the  ejection  period  is  constant. 


C.  Method  I  using  Pressure  P(t)  and  Volume  V ( t ) 

First,  consider  the  case  where  P(t)  and  V (t)  can  be  mea¬ 
sured  directly,  for  example,  with  a  pressure  transducer  and 
a  conductance  catheter,  respectively. 

Define  the  unbiased  ventricular  volume  {v(t))  as 


v(t)  “A1  V(t)  -  Vo  (2) 


Substitution  of  (2)  to  (1)  yields 

m  =  ®  (3) 

v(t) 

Assumption  I  means  that  the  second  derivative  of  E(t) 
is  zero,  that  is, 


E(t) 


2  P(t)v2(t) 
v3(t) 

0 


2P(t)v(t)  _  P(t)v(t)  P(t ) 

v2(t)  v2(t)  v(t) 

(4) 


Solving  (4)  with  respect  to  v(t)  gives  two  solutions.  How¬ 
ever,  the  unbiased  volume  V(t)  —  Vo  must  be  positive  in  the 
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ejection  period,  and  therefore,  the  positive  solution  vp(t)  is  Define  the  ejecting  volume  or  the  integrated  value  {Iff)) 
chosen  as  the  feasible  solution  as  follows:  of  the  outflow  iff)  from  tbe  to  t  as 


vp{t) 


1 

2 m 


2  P(t)v(t)  +  P{t)v{t)— 


(13) 


{2 P{t)v{t)  +  Pft)vff)}2  -  8 Pff)v2ft)Pff) 

(5) 

Rewrite  vpft)  as  vp(P(t),  P(t),  P(t),v(t),v{t)).  From  (2) 
and  (5),  the  dead  volume  Vo  can  be  represented  by 

Vo  =  V(t)~  vp(P(t),  P(t),  P(t),v(t),v(t))  (6) 

In  addition  to  Assumption  I,  suppose  that  Assumption 
II  is  also  established  at  least  during  a  certain  part  of  the 
ejection  period,  and  then 


Since  i(t)  =  0  for  any  t  from  the  end-diastolic  time  {ted)  to 
tbe ,  V(t)  can  be  represented  as 

V{f)  =  V(tbe)  -  I{f)  =  V(ted)  -  I{t)  =  Ved  -  I(t)  (14) 


where  Ved  is  the  end-diastolic  volume.  The  relation  of  (14) 
can  transform  (1)  to 


m 


pjt) 

ved  -V0-  Iff) 


(15) 


Differentiating  V(t)  with  respect  to  t,  we  have 


v(t)  =  Vft) 
v(t)  =  Vft) 

hold.  Thus,  we  have 

Vo  =  V(t)-vp{P{t),P{t),P{t),V(t),V(t)) 


(7)  V{f)  =  - i(t ) 

(8)  Vft)  = -iff) 

Substitution  of  (14),  (16)  and  (17)  to  (10)  yields 


(16) 

(17) 


(9)  Ved-V0(t)d^  I(t)  +  vp{P(t),P(t),P(t),-i{t),-i{t))  (18) 


It  is  true  that  if  both  Assumptions  I  and  II  hold,  the  right 
hand  side  of  (9)  will  be  constant.  However,  the  time  range 
in  which  both  Assumptions  I  and  II  actually  hold  is  not 
exactly  equal  to  the  whole  of  the  ejection  period,  and  then 
the  right  hand  side  of  (9)  must  be  the  function  of  t  in 
the  whole  of  the  ejection  period.  To  avoid  misunderstand, 
the  left  hand  side  of  (9),  Vo  should  be  rewritten  as  a  time 
function  To  it)  as  follows: 


Vj, ft)  d^f  V{f)  -  vp {P{t) ,P(t),  P(t) ,V(t),V(t))  (10) 


An  estimate  (Vo)  of  the  dead  volume  Vo  can  be  obtained 
as  the  value  of  Vo (t)  at  the  time  t  when  ^Vo{t)  =  0,  i.e. , 

Vo  ft)  is  flat  or  a  peak  since  Assumptions  I  and  II  are  si¬ 
multaneously  established  at  that  time.  After  substitution 
of  Vq  to  (1),  the  estimate  ( E{t ))  of  E(t)  can  be  defined  as 


m 


def 


m 

V{t)  -  % 


(ii) 


Finally,  obtain  the  estimate  (E^ax)  of  E^ax  as  follows: 


J-Jm.n.cr.  — 


max  E(t) 

tbe<t<tee 


(12) 


where  tbe  and  tee  are  the  beginning  and  the  end  time  of 
ejection,  respectively. 

D.  Method  II  using  Pressure  P(t)  and  Outflow  i{t ) 

Next,  consider  the  different  case  where  V(t)  cannot  di¬ 
rectly  be  measured  but  ventricular  outflow  (1(f))  can  be 
measured  with  a  flowmeter  or  an  ultrasound  Doppler  im¬ 
age.  Let  P(t)  be  still  measurable. 


In  the  same  way  as  Method  I,  let  an  estimate  ( Ved  —  Vo) 
of  Ved  —  Vo  be  the  value  of  Ved  —  Vo  (t)  at  the  time  t  when 
Ved  —  Vo(t)  is  flat  or  a  peak.  Substitute  Ved  —  Vo  into  (15) 
and  we  have  E(t)  as  follows: 


m 


def 


_P(t) 

Ved -Vo  -m 


(19) 


Finally,  obtain  Emax  according  to  (12). 

E.  Experimental  Protocols  and  Data  Processing 

In  vivo  experiments  were  carried  out  by  using  an  adult 
goat  weighing  50kgs  in  an  open  chest.  An  conductance 
catheter  (Leycom;  Sigma  5)  was  inserted  into  the  left 
ventricle  via  the  apex  to  measure  left  ventricular  volume 
V(t).  Aortic  flow  iff)  was  measured  at  the  ascending  aorta 
with  an  electromagnetic  flow  meter  (Nihon  Kohden;  MFV- 
3100),  and  left  ventricular  pressure  P{t)  was  measured  with 
a  catheter-tip  pressure  transducer  (Camino;  420) . 

In  our  calculation,  a  mathematical  processing  language 
Mathematica  (Wolfram  Research  Inc.)  was  employed  for 
sampled  data  of  P(t),  V(t)  and  i(t)  with  the  sampling  in¬ 
terval  of  At  =  10ms.  In  particular,  its  curve  fitting  func¬ 
tion  Fit  was  used  to  yield  polynomial  functions  P(f),  V(t) 
and  i{t)  in  the  ejection  period  for  eliminating  measurement 
noises.  This  is  because  the  first  and  the  second  derivatives 
of  P(t),  V(t)  and  i(t),  which  are  very  sensitive  to  noise,  are 
required  in  calculation  of  (10)  and  (18). 

III.  Results 

Fig.  1  shows  an  example  of  pressure-volume  loops  (PV- 
loops)  obtained  in  our  experiment.  The  PV-loops  are  de¬ 
picted  during  successive  five  beats  after  manually  clamping 


LW:  V(t)  (mL)  Fig.  4.  Sampled  values  (dots)  of  aortic  flow  rate  lit)  corresponding 

.  .  „  ii-  ,  to  Fig.  lb)  and  its  approximation  (solid  curve), 

big.  1.  An  example  oi  pressure- volume  loops  m  manually  clamping  c 

the  pulmonary  artery. 


the  pulmonary  artery.  The  line  a)  is  an  ESPVR  manually 
determined  on  the  basis  of  the  end-systolic  points.  In  this 
example,  the  corresponding  dead  volume  Vo  can  be  ob¬ 
tained  as  Vo  =  127mL. 

The  sampled  values  of  P(t),  V(t)  and  i(t)  corresponding 
to  the  ejection  period  of  the  loop  b)  in  Fig.  1  are  depicted 
as  dots  in  Figs.  2,  3  and  4,  respectively.  The  sampled  val¬ 
ues  were  approximated  by  the  solid  curves  shown  in  these 
figures.  The  orders  of  the  polynomials  were  3,  3  and  6  for 
P(t),  V(t)  and  i(t),  respectively. 

P(t)  (mmHg) 


Fig.  2.  Sampled  values  (dots)  of  left  ventricular  pressure  P(t)  corre¬ 
sponding  to  Fig.  lb)  and  its  approximation  (solid  curve). 


V(t)  (mL) 


Fig.  3.  Sampled  values  (dots)  of  left  ventricular  volume  V(t)  corr 
sponding  to  Fig.  lb)  and  its  approximation  (solid  curve). 

Fig.  5  shows  Vo(t)  and  Ved  —  Vo(t).  The  upper  curve  Vo(j  , 
was  calculated  from  (10)  substituted  by  approximated  time 
functions  P(t),  V(t)  and  their  derivatives  corresponding  to 
Figs.  2  and  3.  In  the  same  way,  the  lower  curve  Ve([  —  Vo(t) 


was  calculated  from  (18)  substituted  by  P(t),  i(t)  and  their 
derivatives  corresponding  to  Figs.  2  and  4.  In  Fig.  5,  the 
parts  a)  and  b)  of  Vo(t)  satisfy  ifVo(t)  =  0.  However,  the 
part  b)  corresponds  to  the  peak  time  tmax  that  maximizes 
E(t).  Then,  the  estimated  dead  volume  Vo  can  be  given 
as  the  value  of  Vo(t)  of  the  part  a).  In  this  case,  we  have 
Vo  =  126mL.  In  the  same  way,  the  part  d)  of  Ve(i  —  Vo (t) 
can  be  chosen  as  the  estimated  value  Ve(i  —  Vo,  because 
the  part  cl)  is  the  closest  to  the  part  a).  Thus,  we  have 
Ved  -  V0  =  60mL.  ^  _ 

Substitutions  of  Vo  and  Ved  —  Vo  to  (11)  and  (15)  yielded 
the  estimates  E(t)  as  shown  in  Figs.  6  and  7,  respec¬ 
tively.  In  both  figures,  the  shape  of  E(t)  is  nearly  lin¬ 
ear.  The  maximum  value  of  each  E(t)  gives  the  estimated 
maximum  elastance  Emax,  and  thus,  we  have  Emax  = 
2.62mmHg/mL  using  pressure  and  volume  shown  in  Fig.  6, 
and  Emax  =  2.76mmHg/mL  using  pressure  and  flow  rate 
shown  in  Fig.  7. 

The  above  process  was  applied  to  each  ejection  period 
of  the  five  PV-loops  shown  in  Fig.  1.  Table  1  shows  mean 
values  and  standard  deviations  of  estimates  Vo  and  Emax 
averaged  over  the  five  beats.  The  ESPVRs  can  be  depicted 
as  shown  in  Fig.  8  using  the  estimates  of  a)  and  b)  of  Ta¬ 
ble  1.  The  Table  1  and  Fig.  8  indicate  that  the  proposed 
methods  could  yield  the  estimated  values  roughly  close  to 
the  true  values. 


^0  (t)  (mL)  ,  Ved -Vo  (t)  (mL) 


0  0.05  0.1  0.15  0.2 


T ime  ( s ) 

Fig.  5.  Vo  CO  and  Ved  —  Vo(t)  calculated  from  (10)  and  (18),  re¬ 

spectively.  The  peak  values  a)  and  d)  give  the  estimates,  Vq  and 
Ved  —  Vo,  respectively. 


fi(t)=P(t)/{V(t)-'<to} 


Fig.  6.  Estimate  E(t)  of  the  left  ventricular  elastance  E(t)  cal¬ 
culated  from  (11)  on  the  basis  of  approximated  P(t)  and  V{t) 
(solid  curve),  and  their  measured  values  (dots). 


(t)  =  p  (tj/ivecFVo -i  (t) } 


Time  (s) 

Fig.  7.  Estimate  E(t)  of  the  left  ventricular  elastance  E{t)  calculated 
from  (19)  on  the  basis  of  approximated  P(t)  and  i(t)  (solid  curve), 
and  their  measured  values  (dots). 

IV.  Discussion 

The  result  of  Table  1  and  comparison  between  Fig.  1  and 
Fig.  8  imply  that  the  proposed  method  (Methods  I  and  II) 
can  estimate  Emax  beat  by  beat  with  reasonable  accuracy 
without  any  change  in  the  cardiac  load.  However,  the  stan¬ 
dard  deviation  of  each  estimate  is  not  so  small.  This  is 
because  the  method  depends  strongly  on  the  waveform  of 
the  time  series  data  and  is  sensitive  to  noise.  To  cope  with 
this  defect,  the  time  series  data  such  as  P(t),  V(t)  and  i(t) 
should  be  averaged  over  several  beats  which  have  almost 
the  same  cardiac  cycle  as  one  another  in  the  steady  state 
or  under  a  constant  cardiac  load. 

The  proposed  method  does  not  need  the  procedure  for 
optimizing  unknown  parameters  as  the  POM  does,  and 
then  the  proposed  method  requires  much  less  computa¬ 
tional  time  for  estimation  than  the  POM.  Moreover,  the 
POM  must  previously  assign  the  approximation  range 
where  the  elastance  E(t)  can  be  regarded  as  a  line.  On  the 
other  hand,  in  the  proposed  method,  the  adequate  part  in 
which  both  Assumptions  I  and  II  hold  can  automatically 
be  given  as  the  time  when  Vo(f)  or  Ve d  —  Vo  (t)  is  flat  or  a 
peak. 

Senzaki’s  method  [3]  also  does  not  need  the  optimizing 
procedure  but  needs  the  normalized  ventricular  elastance 
curve  2?jv  (ijv)  and  the  recursive  procedure  for  searching 
for  the  peak  time  trnax.  In  this  paper,  unfortunately,  the 
estimation  accuracy  has  not  yet  been  compared  among  the 
proposed  method,  the  POM  and  Senzaki’s  method. 


LW:  V(t)  (inL) 

Fig.  8.  ESPVRs  depicted  by  means  of  the  estimates  of  a)  and  b)  of 
Table  1. 


Table  1  Estimates  Vrj  and  Emax  averaged  over  the  five  beats 


shown 

in  Fig.  1  or  Fig.  8 

True 

value 

a) 

Estimate  using 
P(t)  and  V{t) 

b) 

Estimate  using 
P(t)  and  i(t) 

Vo 

(mL) 

127 

123(±5.3) 

126(-) 

Emax 

(mmHg/mL) 

2.9 

2.5(±0.54) 

2.8(±0.50) 

V.  Conclusions 

The  present  study  has  derived  a  new  method  for  estimat¬ 
ing  the  maximum  ventricular  elastance  Emax  by  assuming 
that  the  following  two  assumptions  are  simultaneously  es¬ 
tablished  in  a  certain  part  of  the  ejection  period:  1)  the 
ventricular  elastance  E(t)  changes  linearly  with  time;  2) 
the  dead  volume  Vo  does  not  change.  The  method  can  es¬ 
timate  Emax  in  a  single  beat  without  any  change  in  preload 
or  afterload  by  measuring  the  ventricular  pressure  and  ei¬ 
ther  the  ventricular  volume  or  outflow.  The  method  does 
not  need  to  previously  assign  the  approximation  range  of 
E(t)  and  does  not  include  the  recursive  procedure  to  search 
for  the  peak  time  tmax  that  gives  Emax.  In  the  experiment 
in  vivo ,  it  could  be  ascertained  that  the  method  can  give 
the  estimate  of  Emax  roughly  close  to  the  true  value  ob¬ 
tained  from  the  conventional  method. 

In  further  studies,  the  estimation  accuracy  of  the  pro¬ 
posed  method  should  be  compared  with  that  of  other 
single-beat  estimation  methods. 
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